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(54) Proportional pressure assist ventilation controlled by a diaphragm electromyographic signal 



(57) A closed loop system uses (a) the intensity of 
the diaphragm electromyogram (EMG) for a given in- 
spiratory volum ; (b) the inspiratory volume for a given 
EMG intensity ; or (c) a combination of (a) and (b) ; in 
view of controlling the level of gas flow, gas volume or 
gas pressure delivered by a mechanical (lung) ventila- 



tor. The closed loop ventilator system enables for auto- 
matic or manual adjustment of the level of inspiratory 
support in proportion to changes in the neuro-ventilatory 
efficiency such that the neural drive remains stable at a 
desired target level. An alarm can also be used to detect 
changes in neuro-ventilatory efficiency in view of per- 
forming manual adjustments. 
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BACKGROUND OF THE INVENTION 

1. Field of the invention: 

[0001] The present invention relates to a system us- 
ing the intensity of the diaphragm electromyogram 
(EMG) at a given lung volume or the lung volume at a 
given EMG intensity to automatically or manually adjust 
the level of inspiratory support in proportion to changes 
in the neuro-ventilatory efficiency. 
[0002] The present invention also relates to a system 
responsive to the intensity of the diaphragm electromy- 
ogram (EMG) measured immediately before the onset 
of inspiratory flow to automatically or manually control 
and maintain an optimum level of extrinsic positive end 
expiratory pressure (PEEP) applied to a patient, and to 
automatically or manually control a duration from the on- 
set of EMG to onset of respiratory flow. 

2. Brief description of the prior art: 

[0003] Prior art algorithms used to create closed-loop 
ventilator systems are based on variables such as tidal 
volume, respiratory rate, inspiratory flow, end-tidal car- 
bon dioxide levels and/or rate of rise in pressure. How- 
ever, none of these parameters can provide a reliable 
measure of the respiratory neural drive because they 
are affected by changes in neuro-mechanical or neuro- 
ventilatory efficiency. 

[0004] Neuro-ventilatory efficiency is a term used to 
express the amount of neural drive (breathing effort) 
needed to obtain a given tidal lung volume. In brief, neu- 
ral drive is converted into mechanical tension, a process 
which is influenced by the muscle length, temperature, 
electrolyte imbalance, etc. The role of inspiratory flow in 
the link between neural drive and mechanical tension 
has previously been suggested; however the proposed 
influence could not be demonstrated for mean inspira- 
tory flow rates up to 1 .4 liters/second. The mechanical 
tension is then translated into pressure, a process which 
is affected by the shape of the diaphragm dome. Finally 
the pressure expands the alveoli and causes air to flow, 
and the translation of pressure to volume depends on 
the elasto-viscous behaviour of the respiratory system. 
Consequently, there are many factors that may influ- 
ence the tidal volume output obtained for a given in- 
crease in neural drive (inspiratory effort). 
[0005] Evaluation of respiratory drive by measure- 
ments such as the rate of rise in pressure or lung volume 
is not reliable when, for example, the muscle length or 
the respiratory system impedance are affected by 
changes in the neuro-ventilatory efficiency. In a patient, 
airway resistance and elastance can change from one 
minute to another and muscle length is continuously al- 
tered. 



OBJECTS AND SUMMARY OF THE INVENTION 

[0006] An object of the present invention is therefore 
to eliminate the drawbacks of the prior art. 
s [0007] Another object of the present invention is to 
provide a closed loop system using: 

( a) the intensity of the diaphragm electromyogram 
(EMG) for a given inspiratory volume; 

10 

(b) the inspiratory volume for a given EMG intensi- 



( c) a combination of ( a) and (b); 

15 

in view of controlling the level of gas flow, gas volume 
or gas pressure delivered by a mechanical (lung) venti- 
lator; the closed loop ventilator system enables for au- 
tomatic or manual adjustment of the level of inspiratory 

20 support in proportion to changes in the neuro-ventilatory 
efficiency such that the neural drive remains stable at a 
desired target level. An alarm can also be used to detect 
changes in neuroventilatory efficiency in view of per- 
forming manual adjustments. 

25 [0008] Another object of the present invention is to 
provide a closed-loop system responsive to the intensity 
of the diaphragm EMG measured immediately before 
the onset of inspiratory flow to quantify pre-inspiratory 
breathing effort in view of automatically or manually ad- 

30 justing a level of extrinsic positive end expiratory pres- 
sure (PEEP) applied to a patient in proportion to chang- 
es in EMG intensity of pre-inspiratory efforts. In this 
manner, the pre-ventilatory intensity of the diaphragm 
EMG can be maintained at a desired, minimum level 

35 such that the pre-inspiratory neural drive remains stable 
at a desired target minimal level. Determination of the 
duration from the onset of EMG to the onset of respira- 
tory flow is also used for quantitative evaluation of the 
intrinsic PEEP, and to guide adjustment of the trigger 

40 sensitivity of the ventilator systems. 

[0009] Different from pressure and ventilatory related 
indexes, the intensity of the EMG represents the tem- 
poral (mean MU (motor unit) rate coding) and spatial 
(MU recruitment) summation of action potentials and is 

4s obtained at the level of the sarcolemma muscle. The in- 
tensity of the EMG is therefore not affected by changes 
in the muscle's neuro-ventilatory coupling. In the 
present invention, the use of crural diaphragm EMG 
rests on the assumption that neural drive to the crural 

so diaphragm is representative for the total respiratory 
drive. It is also based on the condition that neuromus- 
cular transmission and innervation of the crural dia- 
phragm are normal. For breathing with increased de- 
mand this assumption is well founded. Hence, the inten- 
ds sity of the EMG needed to produce a given inspiratory 
volume should express the efficiency relation between 
neural drive and volume output. 
[0010] The objects, advantages and otherfeatures of 
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the present invention will become more apparent upon 
reading of the following non restrictive description of a 
preferred embodiment thereof, given by way of example 
only with reference to the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0011] In the appended drawings: 

Figure 1 is a schematic representation of a set-up 
of an EMG analysis system; 

Figure 2 is a section of oesophageal catheter on 
which an array of electrodes of the EMG analysis 
system of Figure 1 is mounted; 

Figure 3 illustrates a section of oesophageal cath- 
eter on which a second embodiment of the array of 
electrodes is mounted; 

Figure 4 is a graph showing a set of EMG signals 
of the diaphragm (EMGdi signals) detected by pairs 
of successive electrodes of the array of Figure 2; 

Figure 5 is a flow chart showing a method for con- 
ducting a double subtraction technique of the EM- 
Gdi signals; 

Figure 6 is a graph showing the distribution of cor- 
relation coefficients calculated for determining the 
position of the center of the depolarizing region of 
the diaphragm along the array of electrodes of Fig- 
ure 2; 

Figure 7 is a schematic diagram illustrating in the 
time domain a double subtraction technique for im- 
proving the signal-to-noise ratio and to reduce an 
electrode-position-induced filter effect along the ar- 
ray of electrodes of Figure 2; 

Figure 8a is a graph showing the power density 
spectrum of electrode motion artifacts, the power 
density spectrum of ECG, and the power density 
spectrum of EMGdi signals; 

Figure 8b is a graph showing an example of transfer 
function for a filter to be used for filtering out the 
electrode motion artifacts, ECG, and the 50 or 60 
Hz disturbance from electrical mains; 

Figure 9 is a schematic diagram illustrating in the 
frequency domain stabilization by the double sub- 
traction technique of the centerfrequency upon dis- 
placement of the center of the depolarizing region 
of the diaphragm along the array of electrodes of 
Figure 2; 

Figure 1 0 is a schematic block diagram of a system 



66 779 A1 4 

according to the invention for controlling inspiratory 
assist by means of an EMGdi signal obtained with 
the above mentioned double subtraction technique 
and a measurement of the volume of air breathed 
s by the patient by a commercially available system; 

Figure 1 1 is a schematic block diagram of a system 
according to the invention ( a) capable to determine 
the time delay from the onset of EMG to the onset 

10 of inspiratory flow and ( b) using the level of pre- 
inspiratory effort obtained through the EMGdi signal 
intensity (common noise level subtracted) during a 
predetermined time period immediately preceding 
the onset of inspiratory flow to indicate the presence 

15 of "intrinsic PEEP" and to adjust the level of applied 
"extrinsic PEEP" and/or ventilator trigger sensitivity 
such that the level of pre-inspiratory effort is sup- 
pressed, i.e the EMGdi signal intensity (common 
noise level subtracted) during the above mentioned 

20 predetermined time period is close to zero; 

Figure 12a is an exemplary graph of a patient's in- 
spiratory flow versus time for quiet breathing in 
COPD (Chronic Obstructive Pulmonary Disease); 
25 and 

Figure 12b is an exemplary graph of a patient's 
EMG RMS intensity versus time for quiet breathing 
in COPD. 

30 

DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENT 

[0012] Although the preferred embodiment of the 
35 present invention will be described in relation to a double 
subtracted EMGdi signal, it should be kept in mind that 
the concept of the present invention can be used with 
any respiratory muscle signal. 

[0013] To measure EMG activity of the diaphragm 11 
40 (EMGdi) of a human patient 14, an array of electrodes 
such as 12 (Figures 1 and 2) are mounted on the free 
end section 15 of an oesophageal catheter 13, with a 
constant inter-electrode distance d (Figure 2) . As shown 
in Figure 1, the catheter 13 is introduced into the pa- 
45 tient's oesophagus through one nostril orthe mouth until 
the array of electrodes 12 is situated at the level of the 
gastroesophageal junction. The diaphragm 11 and/or 
the oesophagus slightly move during breathing of the 
patient 1 4 whereby the array of electrodes 1 2 also slight- 
so |y moves about the diaphragm 11. As will be explained 
in the following description, automatic compensation for 
this displacement is provided for. 
[0014] According to a preferred embodiment, an elec- 
trode 12 is mounted on the free end section 15 of the 
55 catheter 1 3 by winding stainless steel wire (not shown) 
around that catheter 1 3. The wound stainless steel wire 
presents a rough surface smoothed out by solder, which 
in turn is electroplated with nickel, copper and then gold 
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or silver. Of course, it is within the scope of the present 
invention to use other electrode structures. Also, the 
electrodes 12 can possibly be applied to a nasogastric 
feeding tube (not shown) which is routinely introduced 
in intensive-care unit (ICU) patients, 
[0015] Electric wires (not shown) interconnect each 
pair of successive electrodes such as 1 -7 (Figure 2) with 
a respective one of a group of differential amplifiers 1 6. 
Obviously, these electric wires follow the catheter 13 
from the respective electrodes 12 to the corresponding 
amplifiers 16, and are preferably integrated to the cath- 
eter 13. Preferably, the electric wires transmitting the 
EMGdi signals collected by the various pairs 1 -7 of elec- 
trodes 12 are shielded to reduce the influence of exter- 
nal noise, in particular disturbance from the 50 or 60 Hz 
current and voltage of the electrical mains. 
[0016] The group of differential amplifiers 16 amplifies 
(first subtraction step of a so-called double subtraction 
technique) and band-pass filters each EMGdi signal. 
This first subtraction step may also be carried out in the 
personnel computer 1 9 when the amplifiers 1 6 are sin- 
gle-ended or equivalently designed amplifiers (monop- 
olar readings). 

[0017] In the example illustrated in Figures 1 and 2, 
the free end section 15 of the catheter 13 is provided 
with an array of eight electrodes 1 2 defining seven pairs 
1, 2,3, 4, 5, 6 and 7 of successive electrodes ^respec- 
tively collecting seven different EMGdi signals. Although 
it has been found that EMG activity of the diaphragm 
(EMGdi) can be measured accurately with an oesopha- 
geal catheter 13 provided on the free end section 15 
thereof with an array of eight electrodes 12, a different 
number and/or configuration of pairs of electrodes 12 
can be contemplated depending on the patient's anato- 
my and movement of the diaphragm. Also, the pairs 1 -7 
do not need to be pairs of successive electrodes; as an 
example Figure 3 illustrates an array of nine electrodes 
to form seven overlapping pairs of electrodes 1-7. 
[0018] A major problem in recording EMGdi signals is 
to maintain the noise level as low and as constant as 
possible. Since the electric wires transmitting the EMGdi 
signals from the electrodes 12 to the differential ampli- 
fiers 16 act as an antenna, it is crucial, as indicated in 
the foregoing description, to shield these electric wires 
to thereby protect the EMGdi signals from additional ar- 
tifactual noise. Also, the package enclosing the differ- 
ential amplifiers 16 is preferably made as small as pos- 
sible (miniaturized) and is positioned in close proximity 
to the patient to decrease as much as possible the dis- 
tance between the electrodes 12 and the amplifiers 16. 
[0019] The amplified EMGdi signals are sampled by 
a personal computer 19 through respective isolation 
amplifiers of a unit 1 8, to form signal segments of fixed 
duration. Unit 1 8 supplies electric power to the various 
electronic components of the differential and isolation 
amplifiers while ensuring adequate isolation of the pa- 
tient's body from such power supply. The unit 18 also 
incorporates bandpass filters included in the respective 
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EMGdi signal channels to eliminate the effects of alias- 
ing. The successive EMGdi signal segments are then 
digitally processed into the personal computer 19 after 
analog-to-digital conversion thereof. Th is analog-to-dig- 
s ital conversion is conveniently carried out by an analog- 
to-digital converter implemented in the personal compu- 
ter 19. The personal computer 1 9 includes a monitor 40 
and a keyboard 31 . 

[0020] It is believed to be within the capacity of those 
10 of ordinary skill in the arttoconstructsuitable differential 
amplifiers 16 and an adequate isolation amplifiers and 
power supply unit 18. Accordingly, the amplifiers 16 and 
the unit 1 8 will not be further described in the present 
specification. 

15 [0021] An example of the seven EMGdi signals col- 
lected by the pairs 1 -7 of successive electrodes 1 2 (Fig- 
ures 1 and 2) and supplied to the computer 19 is illus- 
trated in Figure 4. 

[0022] As the diaphragm is generally perpendicularto 

20 the longitudinal axis of the oesophageal catheter 13 
equipped with an array of electrodes 12, only a portion 
of the electrodes 12 are situated in the vicinity of the 
diaphragm. It is therefore important to determine the po- 
sition of the diaphragm with respect to the oesophageal 

25 electrode array. 

[0023] The portion of the crural diaphragm 11 which 
forms the muscular tunnel through which the oesopha- 
geal catheter 1 3 is passed is referred to the "diaphragm 
depolarizing region" (DDR). The thickness of the DDR 

30 is 20-30 mm. It can be assumed that, within the DDR, 
the distribution of active muscle fibers has a center from 
which the majority of the EMGdi signals originate, i.e. 
the "diaphragm depolarizing region center" (DDR cent- 
er). Therefore, EMGdi signals detected on opposite 

35 sides of the DDR center will be reversed in polarity with 
no phase shift; in other words, EMGdi signals obtained 
along the electrode array are reversing in polarity at the 
DDR center. 

[0024] Moving centrally from the boundaries of the 
40 DDR, EMGdi power spectrums progressively attenuate 
and enhance in frequency. Reversal of signal polarity on 
either side of the electrode pair 4 with the most attenu- 
ated power spectrum confirms the position from which 
the EMGdi signals originate, the DDR center. 
45 [0025] Referringto Figure 5, the first task of the com- 
puter 1 9 is to determine the position of the center of the 
DDR along the array of electrodes 12. The center of the 
DDR is repeatedly determined at predetermined time in- 
tervals. 

so [0026] For that purpose, filtering step 505 removes 
from each EMGdi signal the motion artifacts, the elec- 
trocardiogram (ECG) component, and the disturbance 
from the electrical mains. Motion artifacts are induced 
by motion of the electrodes 1 2. More generally, motion 

55 artifacts are defined as a low frequency fluctuation of 
the EMGdi signals' DC level induced by mechanical al- 
terations of the electrode metal to electrolyte interface 
i.e. changes in electrode contact area and/or changes 
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in pressure that the tissue exerts on the electrode. 
[0027] In step 501 , the filtered EMGdi signals from 
step 505 are cross-correlated in pairs. As well known to 
those of ordinary skill in the art, cross-correlation is a 
statistical determination of the phase relationship be- 
tween two signals and essentially calculates the similar- 
ity between two signals in terms of a correlation coeffi- 
cient r (step 502). A negative correlation coefficient r in- 
dicates that the cross-correlated signals are of opposite 
polarities. 

[0028] Figure 6 shows curves of the value of the cor- 
relation coefficient r versus the midpoint between the 
pairs of electrodes from which the correlated EMGdi sig- 
nals originate. In this example, the inter-electrode dis- 
tance is 10 mm. Curves are drawn for distances be- 
tween the correlated pairs of electrodes 12 of 5 mm 
(curve 20), 10 mm (curve 21), 15 mm (curve 22) and 20 
mm (curve 23). One can appreciate from Figure 5 that 
negative correlation coefficients r are obtained when 
EMGdi signals from respective electrode pairs situated 
on opposite sides of the electrode pair 4 are cross-cor- 
related. It therefore appears that the change in polarity 
occurs in the region of electrode pair 4, which is con- 
firmed by the curves of Figure 4. Accordingly, it can be 
assumed that the center of the DDR is situated substan- 
tially midway between the electrodes 12 forming pair 4. 
[0029] For example, the center of the DDR can be pre- 
cisely determined by interpolation (step 503 of Figure 5) 
using a square law based fit of the three most negative 
correlation coefficients of curve 21 obtained by succes- 
sive cross-correlation of the EMGdi signal segments 
from each electrode pair to the EMGdi signal segments 
from the second next electrode pair. Association of the 
center of the DDR to a pair of electrodes 1 2 provides a 
"reference position" from which to obtain EMGdi signal 
segments within the DDR. Such control is essential in 
overcoming the artifactual influence of perpendicular bi- 
polar electrode filtering on the EMGdi power spectrum. 
[0030] It has been experimentally demonstrated that 
EMGdi signals recorded in the oesophagus are satis- 
factory as long as they are obtained from electrode pairs 
(with an inter-electrode distance situated between 5 and 
20 mm) positioned at a distance situated between 5 and 
30 mm on the opposite sides of the DDR center (the 
inter-pair distance being therefore situated between 5 
and 30 mm). Although EMGdi signals obtained from 
these positions offers a clear improvement in accept- 
ance rate, the signal-to-noise ratio during quiet breath- 
ing still tends to remain unsatisfactorily low. The EMGdi 
signal obtained from one electrode pair (for example 
channel 0 in Figure 7) situated in between the two elec- 
trode pairs used to produce the double subtracted sig- 
nal, can be added to this double subtracted signal either 
before as a raw signal or after when RMS or equivalent 
EMGdi signal measure has been computed, in order to 
minimize loss of signal. 

[0031] For example, in Figure 4, the EMGdi signals 
originating from the electrode pairs 3 and 5 situated re- 



spectively 1 0 mm below and 1 0 mm above the DDR are 
strongly inversely correlated at zero time delay. In con- 
trast tothe inversely correlated EMGdi signals, the noise 
components for electrode pairs 3 and 5 are likely to be 
s positively correlated. Hence, as illustrated in Figure 7, 
subtraction of the EMGdi signals 24 and 25 from elec- 
trode pairs 3 and 5 will result into an addition of the cor- 
responding EMGdi signals (signal 26 of Figure 6) and 
into a subtraction, that is an elimination of the common 

io noise components. This technique will be referred to as 
"the double subtraction technique" (step 504 of Figure 
5). Again, the EMGdi signal obtained from one electrode 
pair (for example channel 0 in Figure 7) situated in be- 
tween the two electrode pairs used to produce the dou- 

'5 ble subtracted signal, can be added to this double sub- 
tracted signal either before as a raw signal or after when 
RMS or equivalent EMGdi signal measure has been 
computed, in order to minimize loss of signal. 
[0032] Subtraction step 504 (second subtraction step 

20 of the double subtraction technique) can be carried out 
either in the time domain, or after conversion of signals 
24 and 25 in the frequency domain. Double subtraction 
technique can be performed by subtracting other com- 
binations of signals, for example by subtracting the EM- 

25 Gdi signal segments from electrode pair 2 from the EM- 
Gdi signal segments from electrode pair5 (Figure 4), by 
subtracting signal segments from electrode pair 6 from 
the signal segments from electrode pair 3 and by adding 
these differences, etc. What is important is to subtract 

30 two signals of opposite polarities obtained in the vicinity 
of the muscle. More than two signal pairs of opposite 
polarities can be used in the double subtraction. Again, 
the EMGdi signal obtained from one electrode pair (for 
example channel 0 in Figure 7) situated in between the 

35 two electrode pairs used to produce the double subtract- 
ed signal, can be added to this double subtracted signal 
either before as a raw signal or after when RMS or equiv- 
alent EMGdi signal measure has been computed, in or- 
der to minimize loss of signal. 

40 [0033] The double subtraction technique is carried out 
in step 504 on the pair of EMGdi signals (for example 
the signals from electrode pairs 3 and 5 shown in Figure 
4) identified in step 503, after appropriate filtering of 
these EMGdi signals in step 505. Still again, the EMGdi 

45 signal obtained from one electrode pair (for example 
channel 0 in Figure 7) situated in between the two elec- 
trode pairs used to produce the double subtracted sig- 
nal, can be added to this double subtracted signal either 
before as a raw signal or after when RMS or equivalent 

50 EMGdi signal measure has been computed, in order to 
minimize loss of signal. 

[0034] The graph of Figure 8a shows the power den- 
sity spectrum of the above defined electrode motion ar- 
tifacts, the power density spectrum of ECG, and the 
55 power density spectrum of EMGdi signals. The graph of 
Figure 8b shows an example of transfer function for a 
filter (the dashed line showing the optimal transfer func- 
tion, and the solid linethetransferfunction implemented 
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by the inventors) to be used in step 505 for filtering out 
the electrode motion artifacts, ECG : and the 50 or 60 
Hz disturbance from the electrical mains. Processing of 
the EMGdi signals by the computer 1 9 to follow as close- 
ly as possible the optimal transferfunction of Figure 8b 
will conduct adequately filtering step 505. 
[0035] Therefore, double-subtracted signal segments 
509 are obtained atthe output of step 504 by subtracting 
the EMGdi signal segments from the pair of electrodes 
12 in optimal location above the diaphragm from the EM- 
Gdi signal segments from the pair of electrodes 12 in 
optimal location below the diaphragm. More than two 
signal pairs of opposite polarities can be used in the dou- 
ble subtraction. Again, the EMGdi signal obtained from 
one electrode pair (for example channel 0 in Figure 7) 
situated in between the two electrode pairs used to pro- 
duce the double subtracted signal, can be added to this 
double subtracted signal either before as a raw signal 
or after when RMS or equivalent EMGdi signal measure 
has been computed, in order to minimize loss of signal. 
[0036] Referring back to Figure 5, step 506 calculates 
the RMS (root-mean-square) or equivalent or similar 
value 51 0 of the double-subtracted signal segments 509 
produced in step 504. The increase in intensity obtained 
with the double subtraction technique is associated with 
a twofold increase in RMS values. RMS values obtained 
with the double subtraction technique are closely and 
linearly related to the original signals. It should be kept 
In mind that the RMS value can be replaced by any other 
value representative of the strength of the double-sub- 
tracted signal segments 509. 

[0037] The digital RMS signal segment value 510 cal- 
culated by the computer 1 9 in step 506 is finally digital- 
to-analog converted to an on-line analog RMS value 508 
(step 507) in view of controlling a lung ventilator 54 (Fig- 
ure 1 0). It should be mentioned that it is within the scope 
of the present invention to supply a digital value 508. 
[0038] The double subtraction technique compen- 
sates for the changes in signal strength and frequency 
caused by movement of the diaphragm 11 (Figure 1) 
and/or the oesophagus during breathing of the patient 
14 causing movement of the array of electrodes 12 with 
respect to the diaphragm 11 . Referring to Figure 9, off 
center of the array of electrodes 12 (electrode-position- 
induced filter effect) causes a variation of center fre- 
quency values due to filtering (see curves 27 and 28) 
for the EMGdi signals from the electrode pairs 3 and 5. 
The double subtraction technique eliminates such vari- 
ation of center frequency values as indicated by curve 
29 as well as variation of signal strength. Therefore, the 
reciprocal influence of the position of the DDR center on 
the EMGdi signal frequency content is eliminated by the 
double subtraction technique. 

[0039] It has been found that the double subtraction 
technique may improve the signal-to-noise ratio by more 
than 2 dB and reduce an electrode-position-induced fil- 
ter effect. Double subtraction technique is also respon- 
sible for a relative increase in acceptance rate by more 



than 30%. 

[0040] Noise of non diaphragmatic origin or artifactual 
signals are strongly correlated at zero time delay and 
equal in polarity between all pairs of electrodes 12. 

s Hence, this noise of non diaphragmatic origin or artifac- 
tual signals appear as a common mode signal for all 
electrode pairs and therefore, are substantially reduced 
by the double subtraction technique. 
[0041] In the following description, it should be con- 

10 sidered that the flow and volume of air breathed by the 
patient can be measured by any commercially available 
system. 

Neuro-ventllatory efficiency: 

15 " ' ~ 

[0042] The neuro-ventilatory efficiency is obtained by 
relating the diaphragm EMGdi signal intensity to chang- 
es in lung volume, or by relating the lung volume to 
changes in diaphragm EMGdi signal intensity. Since the 
20 relationship between the diaphragm EMGdi signal in- 
tensity and the lung volume is not linear, this non-line- 
arity is minimized by expressing: 

the intensity of the diaphragm EMGdi signal for a 
25 given volume change from end-expiratory lung vol- 
ume, for example the EMGdi signal intensity ob- 
tained during 400 ml inspiration starting from end- 
expiratory lung volume (in the present disclosure, 
Intensity is intended to encompass the mean, peak, 
30 median and total RMS intensity of the diaphragm 
EMGdi signal); or 

the lung volume obtained at a given diaphragm EM- 
Gdi signal intensity. 

35 

A relatively small tidal lung volume is suitable because 
the relationship between diaphragm EMGdi signal in- 
tensity and lung volume is relatively linear at this low 
range. Secondly, the use of a fixed, given tidal volume 
40 or diaphragm EMGdi signal intensity will protect against 
the non-linear influences and allows for a reliable esti- 
mation of relative changes in neuro-ventilatory efficien- 
cy. 

[0043] In this manner, a ventilatory efficiency index 
45 expressing: 

the EMGdi signal intensity for a given inspiratory 
lung volume starting from the end-expiratory lung 
volume; or 

50 

the lung volume for a given diaphragm EMGdi sig- 
nal intensity; 

is calculated. If the EMGdi signal intensity for the above 
55 mentioned given inspiratory lung volume or the lung vol- 
ume for the above mentioned given diaphragm EMGdi 
signal intensity is changing, the above indicated index 
will also change and this change can be expressed in 



11 



EP 1 366 779 A1 



12 



percentage (%). For example, using the diaphragm EM- 
Gdi signal intensity for the above mentioned fixed, given 
inspiratory lung volume, an increased EMGdi signal in- 
tensity for the above mentioned given inspiratory lung 
volume will increase the index but will express a reduc- 
tion in the neuro-ventilatory efficiency, and a decreased 
EMGdisignal intensity forthat given inspiratory lung vol- 
ume will reduce the index but will express an improve- 
ment of the neuro-ventilatory efficiency. 
[0044] In the following description, an example using 
the EMGdi signal intensity for a fixed, given inspiratory 
lung volume will be given. However, it is within the scope 
of the present invention to use the lung volume for a 
fixed, given diaphragm EMGdi signal intensity. 
[0045] Referring now to Figure 1 0 a preferred, practi- 
cal embodiment is described. A neuro-ventilatory effi- 
ciency computation device 601 receives the signal 508 
of Figure 5 as well as the given, fixed inspiratory lung 
volume. Device 601 comprises a unit 602 for determin- 
ing the intensity of the signal 508 for the given inspiratory 
lung volume. Although it is not illustrated, it is within the 
scope of the present invention to calculate, in unit 602, 
the peak, mean , median or any other intensity measure 
of signal 508 for the given inspiratory lung volume. If the 
intensity of signal 508 for the given inspiratory lung vol- 
ume has increased at least by a given percentage (step 
603), i.e. the neuro-ventilatory efficiency index has in- 
creased at least by said given percentage, the pressure, 
flow, or volume assist unit 604 is controlled by a unit 606 
in view of increasing the magnitude of the pressure as- 
sist to the patient by a preset increment until the intensity 
of the signal 508 for the given inspiratory lung volume 
is restored to a predetermined, preset value. 
[0046] Still referring to Figure 10, if the intensity for 
the given inspiratory lung volume has decreased at least 
by a given percentage (step 607), i.e. the neuro-venti- 
latory efficiency index has decreased at least by said 
given percentage, the pressure assist unit 604 is con- 
trolled by the unit 608 in view of decreasing the magni- 
tude of the pressure assist by a preset increment until 
the intensity of the signal 508 for the given inspiratory 
lung volume is restored to the predetermined, preset 
value. Although it is not illustrated, it is within the scope 
of the present invention to calculate, in unit 602, the 
peak, mean, median or any other intensity measure of 
signal 508 for the given inspiratory lung volume, instead 
of the intensity of this signal. Also, the signals at the out- 
puts of the units 606 and 608 can be used to generate 
an alarm orto manually adjust the pressure, flow or vol- 
ume assist to the patient. 

[0047] The response time is adjustable. The time 
base used to calculate trends in the EMG intensity for a 
given volume or vice versa and used forthe corrections 
is relatively slow (minutes) and the levels of applied sup- 
port can be limited within a safe range. Again, an alarm 
can be generated or the pressure assist can be manu- 
ally or automatically adjusted. 

[0048] The pressure, flow, or volume assist unit 604 



can be any device which can be controlled to generate 
any airway pressure of adjustable magnitude, for exam- 
ple any source of compressed gas, or a flow or volume 
pump. Of course, airway 605 refers to or, to the least, 

5 includes the patient's respiratory airway. 

[0049] In this manner, the pressure assist unit 604 
provides a pressure, flow, or volume assist that is ad- 
justed in proportion to changed in neuro-ventilatory ef- 
ficiency which is the EMGdi signal intensity at a given 

10 lung volume or vice versa. The pressure, flow, or volume 
assist unit continuoulsy operates to maintain a tracheal 
pressure, flow or volume that is adjusted in proportion 
to changes in neuro-ventilatory efficiency which is the 
EMGdi signal intensity at a given lung volume or vice 

is versa. 

Pre-lnsplratory breathing effort: 

[0050] Acommon problem with mechanically ventilat- 

2Q ed patients is that the patients' inspiratory effort will not 
immediately cause an inspiratory airflow so called "in- 
trinsic PEEP" or "auto PEEP" which leads to a decrease 
in the neuro-ventilatory efficiency. The effect of "intrinsic 
PEEP" can be counteracted by the application of an "ex- 

25 trinsic PEEP" . However, there are no easy applicable 
techniques to determine when the applied level of "ex- 
trinsic PEEP" is adequate. The level of pre-inspiratory 
effort obtained through the EMGdisignal intensity (com- 
mon noise level subtracted) during for example a 100 

30 milliseconds (ms) period immediately preceding the on- 
set of inspiratory flow can be used to indicate the pres- 
ence of "intrinsic PEEP", and the level of applied "ex- 
trinsic PEEP" can be adjusted such thatthe level of pre- 
inspiratory effort is suppressed i.e the EMGdi signal in- 

35 tensity (common noise level subtracted) during the 
above mentioned 1 00 ms period before onset of inspir- 
atory flow is close to zero. A feedback loop can then be 
used to maintain the level of pre-inspiratory effort sup- 
pressed by adjusting as explained above the level of 

40 "extrinsic PEEP". 

[0051] Just a word to mention that the above men- 
tioned period of 1 00 ms can be replaced by a longer or 
shorter time period immediately preceding the onset of 
inspiratory flow or by the neuro-ventilatory delay 800 

■*s (Figure 12b), i.e. the time period between the onset of 
EMG 801 (Figure 1 2b) and the onset of inspiratory flow 
802 (Figure 12a). 

[0052] Figure 11 of the appended drawings illustrates 
a preferred, practical embodiment 700. 
50 [0053] In the embodiment 700, an integrator 71 3 is re- 
sponsive to the RMS EMG signal 508 to continuously 
calculate the EMG intensity for the above mentioned 
100 ms period or neuro-ventilatory delay 800. 
[0054] Embodiment 700 also comprises an inspirato- 
rs ry flow detector 702 responsive to the patient's inspira- 
tory flow 703 measured, as indicated in the foregoing 
description, through any commercially available sys- 
tem, to produce an output signal 705 representative of 
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EMG activity. 

[0055] The embodiment 700 of Figure 11 also com- 
prises a neuro-ventilatory delay calculator 704 respon- 
sive to ( a) the detection of a RMS EMG signal intensity 
higher than the common noise level (5%). and (b) the 
detection of the onset of inspiratory flow by the detector 
702 to calculate the neuro-ventilatory delay 800 (Figure 
12b). 

[0056] A detector 71 4 is responsive to the EMG inten- 
sity calculated by the integrator 713 to detect the level 
of EMG intensity 803 (Figure 12b) at the onset of inspir- 
atory flow 802 (Figure 1 2a) to trigger an alarm 71 6 when 
the level of the EMG intensity 803 at the onset of inspir- 
atory flow 802 is higherthan a given limit (detector 71 5). 
Upon triggering of the alarm 716, the level of applied 
"extrinsic PEEP" is either automatically or manually in- 
creased (device 708). 

[0057] The detector 714 is responsive to the EMG in- 
tensity calculated by the integrator 71 3 to detect the lev- 
el of EMG intensity 803 (Figure 12b) at the onset of in- 
spiratory flow 802 (Figure 12a) to trigger an alarm 720 
when the level of the EMG intensity 803 at the onset of 
inspiratory flow 802 is lower than a given limit (detector 
719). Upon triggering of the alarm 720, the level of ap- 
plied "extrinsic PEEP" is either automatically or manu- 
ally decreased (device 711). 

[0058] It should be mentioned that feedback from the 
neuro-ventilatory delay or pre-inspiratory EMG activity 
can also be used to adjust the sensitivity of the ventila- 
tors trigger functions. 

[0059] Again, the time base used for these corrections 
is preferably relatively slow (minutes) and the levels of 
"extrinsic PEEP" can be limited within a safe range. 
[0060] The pressure assist unit 604 can be any device 
which can be controlled to generate any airway flow and/ 
or pressure of adjustable magnitude, for example any 
source of compressed gas, or a flow or volume pump. 
[0061] In this manner, the delay from the beginning of 
the mechanically ventilated patients' inspiratory effort to 
the onset of the inspiratory assist will be minimized. 
[0062] Although the present invention has been de- 
scribed hereinabove with reference to preferred embod- 
iments thereof, these embodiments can be modified at 
will, within the scope of the appended claims, without 
departing from the spirit and nature of the subject inven- 
tion. 



Claims 

1 . A neuro-ventilatory efficiency computation method 
for monitoring/controlling the level of ventilatory as- 
sist produced by a ventilatory assistance system, 
comprising: 

a) receiving a first signal representative of in- 
spiratory effort and having a first amplitude; 

b) receiving a second signal representative of 



a lung volume and having a second amplitude; 
c)calculating a relation between said first and 
second amplitudes at predetermined intervals; 
and 

s d) increasing or decreasing the ventilatory as- 

sist level depending on whether a present cal- 
culated value of said relation is higher or lower 
than a past calculated value of said relation by 
an amount exceeding a given threshold. 

10 

2. A neuro-ventilatory efficiency computation method 
as in claim 1 , wherein said relation calculating com- 
prises calculating a ratio between said first and sec- 
ond amplitudes at predetermined time intervals. 

3. A neuro-ventilatory efficiency computation method 
as in claim 1 , wherein said relation calculating com- 
prises calculating a ratio between said first and sec- 
ond amplitudes at intervals when one of said first 

20 and second amplitudes reaches a predetermined 
level. 

4. A neuro-ventilatory efficiency computation method 
as in claim 1 , 2 or 3, wherein the ventilatory assist 

25 level increasing or decreasing comprises increas- 
ing the ventilatory assist level when said present 
calculated value of said relation is higher than said 
past calculated value of said relation by an amount 
exceeding the given threshold, and decreasing the 

30 ventilatory assist level when said present calculated 
value of said relation is lower than said past calcu- 
lated value of said relation by an amount exceeding 
said given threshold. 

35 5. A neuro-ventilatory efficiency computation method 
as in claim 1,2.3 or4, wherein receiving the second 
signal representative of a lung volume comprises 
receiving a signal representative of a given lung vol- 
ume. 

6. A neuro-ventilatory efficiency computation method 
as in claim 1 , 2, 3, 4 or 5, wherein receiving the first 
signal representative of inspiratory effort comprises 
receiving a signal representative of a given level of 

*5 inspiratory effort. 

7. A neuro-ventilatory efficiency computation method 
as in claim 1 , 2, 3, 4, 5 or 6, further comprising gen- 
erating an alarm signal when said present calculat- 

so ed value of said relation is higher or lower than the 
past calculated value of said relation by an amount 
exceeding said given threshold. 

8. A neuro-ventilatory efficiency computation method 
55 as in claim 1 , 2, 3, 4, 5, 6 or 7, comprising manually 

performing said increasing or decreasing of the 
ventilatory assist level. 
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9. A neuro-ventilatory efficiency computation method 
as in claim 1, 2, 3, 4, 5, 6, 7 or 8, comprising ex- 
pressing the first signal representative of inspiratory 
effort as one of the following values: a mean of said 
first amplitude, a median of said first amplitude, and 
a peak of said first amplitude. 

10. A neuro-ventilatory efficiency computation method 
as in claim 1, 2, 3, 4, 5, 6, 7 or 8, comprising ex- 
pressing the second signal representative of a lung 
volume as one of the following values: a mean of 
said second amplitude, a median of said second 
amplitude, and a peak of said second amplitude. 

11. A neuro-ventilatory efficiency computation method 
as in claim 1, 2, 3, 4, 5, 6, 7, 8, 9 or 10, wherein 
receiving the first signal representative of inspirato- 
ry effort comprises receiving an electromyographic 
signal from at least one muscle of a patient. 

12. A neuro-ventilatory efficiency computation device 
for monitoring/controlling the level of ventilatory as- 
sist produced by a ventilatory assistance system, 
comprising: 

a) a first input for receiving a first signal repre- 
sentative of inspiratory effort and having a first 
amplitude; 

b) a second input for receiving a second signal 
representative of a lung volume and having a 
second amplitude; 

c) means for calculating a relation between said 
first and second amplitudes at predetermined 
intervals; and 

d) means for increasing or decreasing the ven- 
tilatory assist level depending on whether a 
present calculated value of said relation is high- 
er or lower than a past calculated value of said 
relation by an amount exceeding a given 
threshold. 

13. A neuro-ventilatory efficiency computation device 
as in claim 12, wherein: 

the calculating means comprises a divider re- 
sponsive to the first and second amplitudes for 
calculating a ratio between said first and sec- 
ond amplitudes at predetermined intervals; 
the increasing or decreasing means comprises: 

a comparator responsive to the present 
calculated value and the past calculated 
value of said relation for producing a signal 
representative of a comparison between a 
present calculated value of said relation 
and a past calculated value of said relation; 
an adder interposed between the compa- 
rator and the ventilatory assistance system 



for adding a preset increment to or sub- 
tracting a preset decrement from said ven- 
tilatory assist level when the comparison 
signal exceeds a given threshold. 

5 

14. A neuro-ventilatory efficiency computation device 
as in claim 1 2 or 1 3, wherein said calculating means 
comprises means for calculating said relation at 
predetermined time intervals. 

10 

15. A neuro-ventilatory efficiency computation device 
as in claim 1 2 or 1 3, wherein said calculating means 
comprises means for calculating said relation at in- 
tervals when one of said first and second ampli- 

'5 tudes reach a predetermined level. 

16. A neuro-ventilatory efficiency computation device 
as in any one of claims 1 3 to 1 5, wherein said adder 
comprises means for adding said preset increment 

20 to said ventilatory assist level when said present 
calculated value of said relation is higher than said 
past calculated value of said relation by an amount 
exceeding said given threshold, and means for sub- 
tracting said preset decrement from said ventilatory 

25 assist level when said present calculated value of 
said relation is lowerthan said past calculated value 
of said relation by an amount exceeding said given 
threshold. 

30 17. A neuro-ventilatory efficiency computation device 
as in any one of claims 1 2 to 1 6, wherein the second 
signal representative of a lung volume is a signal 
representative of a given lung volume. 

35 18. A neuro-ventilatory efficiency computation device 
as in any one of claims 12 to 1 6, wherein the second 
signal representative of inspiratory effort is a signal 
representative of a given level of inspiratory effort. 

40 19. A neuro-ventilatory efficiency computation device 
as in any one of claims 1 2 to 1 B, further comprising 
an alarm generator to produce an alarm signal 
when said present calculated value of said relation 
is higher or lower than the past calculated value of 

45 said relation by an amount exceeding said given 
threshold. 

20. A neuro-ventilatory efficiency computation device 
as in any one of claims 12 to 1 9, wherein said adder 

so comprises a manual adjustment system to add said 
preset increment to or subtracting said preset dec- 
rement from said ventilatory assist level. 

21. A neuro-ventilatory efficiency computation device 
55 as in anyone of claims 12 to 20, comprising means 

for expressing the first signal representative of in- 
spiratory effort by means of one of the following val- 
ues: a mean of said first amplitude, a median of said 
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first amplitude, and a peak of said first amplitude. 

22. A neuro-ventilatory efficiency computation device 
as in any one of claims 12 to 20, further comprising 
means for expressing the second signal represent- s 
ative of a lung volume by means of one of the fol- 
lowing values: a mean of said second amplitude, a 
median of said second amplitude, and a peak of 
said second amplitude. 

23. A neuro-ventilatory efficiency computation device 
as in any one of claims 12 to 22, wherein the first 
signal representative of inspiratory effort is an elec- 
tromyographic signal from at least one muscle of a 
patient. ?5 
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